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Abstract 

We prove that, contrary to the common behef, the classical Maxwell electrodynam- 
ics of a point-like particle may be formulated as an infinite-dimensional Hamiltonian 
system. We derive well defined quasi-local Hamiltonian which possesses direct physical 
interpretation being equal to the total energy of the composed system (field -|- particle). 
The phase space of this system is endowed with an interesting symplectic structure. 
We prove that this structure is strongly non-degenerated and, therefore, enables one 
to define the consistent Poisson bracket for particle's and field degrees of freedom. We 
strees that this formulation is perfectly gauge-invariant. 
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1 Introduction 



Classical electrodynamics of charged, point-like particles is usually based on the celebrated 
Lorentz-Dirac equation which, although very useful in many applications, has certain 
inherent difficulties. There have been many attempts to solve this problem (see e. g. ^ - 
for the rewiev). However, there is no solution which is accepted by all physicists. 

The aim of this paper is to show that despite of these problems the classical electrodynam- 
ics of point particles have well defined Hamiltonian structure. We would like to stress that 
the existence of this structure is nontrivial (up to our knowledge, it is completely unknown). 

The derivation of the Lorentz-Dirac equation is based on the decomposition of the elec- 
tromagnetic field into "retarded" and "incoming" components: 



In this approach it is impossible to correctly formulate the initial value problem and find 
the generator of time translations (see 0). The dynamics of external "incoming" field is 
already given and the retarded component is fully determined by the particle's motion. 
Therefore, the field degrees of freedom are completely eliminated and there is simply no 
room for the Hamiltonian description. However, recently it was shown (cf. @]) that classical 
electrodynamics of point-like objects may be formulated as an infinite-dimensional dynamical 
system. In this approach particles and field degrees of freedom are kept on the same footing. 
There is no equations of motion for the particle. It was shown that the conservation law 
for the total (renormalized) four-momentum is equivalent to a certain boundary condition 
for the electromagnetic field along the particle's trajectory. Together with this condition, 
the theory (called electrodynamics of moving particles) becomes causal and complete: initial 
data for both the field and the particles uniquely imply the evolution of the system. This 
means e. g. that the particles trajectories may also be calculated uniquely from the initial 
data. Because this approach is relatively little known (and it is crucial for the present paper) 
we present a short review in the next Section. 

It turns out (cf. [^) that the electrodynamics of moving particles possesses very inter- 
esting Lagrangian structure. Let us note that the standard variational principle used in 
electrodynamics cannot be extended to the theory containing also point-like particles inter- 
acting with the electromagnetic field. Such a principle is based on the following Lagrangian, 
written usually in textbooks: 



J ret ' J in 




'total Li j^dxwell ~l~ -^particle ~l~ Lij^t i 




with 





Lparticle ■ Tn6^ , 



(1.4) 



and the interaction term given by 



'int ■ 



(1.5) 
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Here by 5^ we denote the Dirac delta distribution localized on the particle trajectory (. The 
above Lagrangian may be used to derive the trajectories of the test particles, when the field 
is given a priori. In a different context, it may also be used to derive Maxwell equations, 
if the particle trajectories are given a priori. Simultaneous variation with respect to both 
fields and particles leads, however, to a contradiction, since the Lorentz force will be always 
ill defined due to Maxwell equations. 

But already in the context of the inhomogeneous Maxwell theory with given point-like 
sources, the variational principle based on Lagrangian (|1.2|) is of very limited use, since the 
interaction term L^t becomes infinite. As a consequence, the Hamiltonian of such a theory 
will always be ill defined, although the theory displays a perfectly causal behaviour. 

To improve this bad feature of the theory, a new, quasi-local variational principle for 
the Maxwell field with given sources was proposed in 0. This new variational principle is 
based on the quasi- local Lagrangian which, contrary to ( [L.2| ), is well defined, i.e. produces 
no infinities. It was proved (see 0) that adding to that Lagrangian the particle Lagrangian 
( |1.4| ) and varying it with respect to both fields and particles is now possible and does not 
lead to any contradiction. As a result, one obtains precisely the electrodynamics of moving 
particles proposed in p. 

In the present paper we give the Hamiltonian formulation of the above theory, i. e. 
we shall prove that electrodynamics of moving particles may be formulated as an infinite- 
dimensional Hamiltonian system. For the simplicity we consider here only one particle case. 
It is of course possible to generalise this result to many particles, however, it is technically 
much more complicated. 

Obviously, the Hamiltonian description based on the standard Lagrangian (|1.2|) is incon- 
sistent in the case of a point-like sources. In particular, the Hamiltonian of such a theory is 
always ill defined. Our approach, based on the quasi-local Lagrangian defined in [0] leads via 
an appropriate Legendre transformation to the well defined quasi-local Hamiltonian struc- 
ture for fields interacting with a charged particle. In particular, the Hamiltonian of the 
composed "particle + field" system equals numerically to the total energy of the composed 
system and, therefore, has a direct physical interpretation. 

Moreover, it turns out that the above theory possesses highly nontrivial Poisson bracket 
structure. This structure is defined via the symplectic form living in the phase space of 
the entire system (particle + field). The characteristic feature of the infinite-dimensional 
symplectic manifold is that the symplectic 2-form is in general only weakly nondegenerate, 
cf. [|15|. It means that there need not exist a Hamiltonian vector field Xyr corresponding to 
every given functional JF on the phase space of the system. For example, our system has well 
defined Hamiltonian, however, corresponding vector field generated dynamics is not defined 
throughout the phase space. We prove (cf. Theorem |) that this vector field is well defined 
if and only if we restrict the phase space to the states fulfilling the fundamental equation for 
electrodynamics of moving particles. 

Our main result consists in the Theorem ^ which says, that on the reduced phase 
space the symplectic 2-form becomes strongly nondegenerate. This nice mathematical result 
enables one to define the Poisson bracket structure in the space of functionals over the 
reduced phase space. To our knowledge the above Poisson bracket structure is the first fully 
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consistent structure for the theory of interacting particles and fields, i. e. when the particles 
and fields variables are kept at the same footing. 
The paper is organized as follows. 

Section ^ contains the main results of the new approach to classical electrodynamics of 
point particles proposed in ^j. 

In Section § we present a new technique, developed in [0], which enables us to describe 
at the same footing the field and the particle's degrees of freedom. For this purpose we 
formulate the Hamiltonian structure of any relativistic, hyperbolic field theory with respect 
to a non-inertial reference frame defined as a rest-frame for an arbitrarily moving observer. 

In Section § we show how to extend the above approach to the case of electrodynamical 
field interacting with point particles. 

For the reader convenience we present in Section |^ a new Hamiltonian structure for 
Maxwell electrodynamics. The new electrodynamical Hamiltonian corresponds to a sym- 
metric energy-momentum tensor and it is related via a simple Legendre transformation with 
electrodynamical Lagrangian derived in 

Next three Sections present consistent Hamiltonian structure for electrodynamics of a 
point particle, i. e. we derive the quasi-local Hamiltonian and find the well defined formula 
for the Poisson bracket. In Section ^ we show that the above Poisson bracket structure is 
consistent with Poincare algebra structure of relativistic theory. 

Finally, in Section |1^ we present the Hamiltonian formulation for the particle interacting 
not only with the radiation field, but also with a fixed, external potential, produced by a 
heavy external device. This is a straightforward extension of our theory, where the homoge- 
neous boundary condition for the radiation field is replaced by an inhomogeneous condition, 
the inhomogeneity being provided by the external field. 

2 The new approach to electrodynamics of a point par- 
ticle 

In the present Section we briefly sketch the new approach to electrodynamics of a point 
particle presented in |^. 

Let y = q(t) with t = y^, be the coordinate description of a time-like world line ( of a 
charged particle with respect to a laboratory frame, i.e. to a system (y^), /i = 0, 1,2,3; of 
Lorentzian space-time coordinates. 

The theory contains as a main part the standard Maxwell equations with point-like 
sources: 

d[\ffj.u] = , 

d,r = eu-'5c, (2.1) 

where stands for the particle four-velocity and 5^ denotes the ^-distribution concentrated 
on the smooth world line (: 
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Here v = {v'^) = (q^) is the corresponding 3-velocity and denotes the square of its 3- 
dimensional length (we use the Heaviside-Lorentz system of units with c = 1). In the case 
of many particles the total current is a sum of contributions corresponding to many disjoint 
world lines and the value of charge is assigned to each world line separately. 

For a given particle's trajectory, equations ( p.l|) define a deterministic theory. This means 
that initial data for the electromagnetic field uniquely determine its evolution. However, if 
we want to treat also the particle's initial data (q,v) as dynamical variables, the theory 
based on the Maxwell equations alone is no longer deterministic: the particle's trajectory 
can be arbitrarily modified in the future or in the past without changing the initial data. 

This non-completeness of the theory is usually attributed to the fact that the particle's 
equations of motion are still missing. However, it was proved in that the field initial data 
fully determine the particle's acceleration and this is due to Maxwell equations only, without 
postulating any equations of motion. More precisely, there is a one-to-one correspondence 
between the (r~^)-term of the field in the vicinity of the particle and the acceleration of the 
particle. The easiest way to describe this property of Maxwell theory is to use the particle's 
rest-frame. For this purpose consider the 3-dimensional hyperplane orthogonal to ( at 
the point (t, q(t)) G (. We shall call the "rest frame hyperplane". Choose on any 
system (x*) of cartesian coordinates centered at the particle's position and denote by r the 
corresponding radial coordinate. The initial data for the field on are given by the electric 
induction field D = {D^) and the magnetic induction field B = (i?*) fulfilling the conditions 
div B = and div D = e (5o . Maxwell equations can be solved for arbitrary data, fulfilling 
the above constraints, but the solution will be usually non-regular, even far away from the 
particles. To avoid singularities propagating over a light cone from (t, q(t)), the singular 
part of the data in the vicinity of the particle has to be equal to 



where a = {a^) is the acceleration of the particle (in the rest frame we have a° = 0) and 
0(1) denotes the nonsingular part of the field (the magnetic field B^{r) cannot contain any 
singular part). This gives the one-to-one correspondence between the (r~^)-term of the field 
and the particle's acceleration, which is implied by the regularity of the field outside of the 
trajectory 

Hence, for regular solutions, the time derivatives (D,B,q, v) of the Cauchy data 
(D, B, q, v) of the composed (fields + particle) system are uniquely determined by the data 
themselves. Indeed, D and B are given by the Maxwell equations, q = v and v may be 
uniquely calculated from equation ( p.2| ). Nevertheless, the theory is not complete and its 
evolution is not determined by the initial data. This non-completeness may be interpreted 
as follows. Field evolution takes place not in the entire Minkowski space M, but only outside 
the particle, i. e. in a manifold with a non-trivial boundary Mq := M — {(}. The boundary 
conditions for the field are still missing. 

To find this missing condition, the following method was used. Guided by an extended 
particle model, an "already renormalized" formula was proposed in P|, which assigns to each 
point {t,q^{t)) of the trajectory a four-vector p^{t) according to 




(2.2) 
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DEFINITION 1 

px{t) := muxit) + P jjT\ - T^t)\) dT.^ , (2.3) 

where denotes the symmetric energy-momentum tensor of the Maxwell field and T{t)^ 
denotes the energy-momentum tensor corresponding to the electromagnetic field produced by 
the uniformly moving particle along the stright line tangent to the trajectory ( at (t, q(t)). 
"P" denotes the principal value of the singular integral and S is an arbitrary space-like 
hypersurface passing through the point (t,q(t)). 

It was proved in that the four-vector ( |2.3|) is well defined and does not depend on 
the particular choice of S. It is interpreted as the total four-momentum of the physical 
system composed of both the particle and the field. For a generic trajectory ( and a generic 
solution of Maxwell equations ( ^TTP this quantity is not conserved, i. e. it depends upon t. 
The conservation law 

|p^W=0 (2.4) 

is proposed as an additional equation, which completes the theory. It was shown that, due to 
Maxwell equations, only 3 among the 4 equations ( p.4|) are independent. Given a laboratory 
reference frame, one may take e. g. the conservation of the momentum p = (p'^): 

|p = (2.5) 

as independent equations. They already imply the energy conservation ^P^{t) = 0. In 
Section ^ we will show, that the above momentum p is equal to the momentum canonically 
conjugate to the position of the particle, whereas po is equal to the total Hamiltonian of the 
composed (particle + field) system. 

It has been proved in that, due to Maxwell equations, the integral (global) condition 
(|2.5| ) is equivalent to a (local) boundary condition for the behaviour of the Maxwell field 
in the vicinity of the trajectory. The condition was called the fundamental equation of the 
electrodynamics of moving particles. In particle's reference frame it may be formulated as 
a relation between the {r~^) and the (r°) terms in the expansion of the radial component of 
the electric field in the vicinity of the particle: 



47r Vr^ r 



)+/3 + 0(r), (2.6) 



where by 0(r) we denote terms vanishing for r —>■ like r or faster. For a given value of 
r both sides of (|2.(j|) are functions of the angles (only the term is angle-independent). 



The relation between the acceleration and the (r ) - term of the electric field given in (|2.2| ) 
may be rewritten in terms of the component a of this expansion: 

a = —CQi — (2.7) 
r 

(it implies that the quadrupole and the higher harmonics of a must vanish for regular 
solutions). One can prove the following theorem (see H): 
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THEOREM 1 The conservation law (2^) is equivalent to the following boundary condition 

DP{ma + e^(3) = , (2.8) 
where DP(f) denotes the dipole part of the function f on the sphere S'^. 
The main result of ||^ may be summarised in the following 

THEOREM 2 Maxwell equations together with the boundary condition ( ^.8^ ) define com- 
plete, causal and fully deterministic theory: initial data for particles and fields uniquely 
determine the entire history of the system. 



3 Hamiltonian description of a relativistic field theory 
in the co-moving frame 

3.1 Definition of the co- moving frame 

To construct the Hamiltonian formulation of the above theory, we will need a Hamiltonian 
description of electrodynamics with respect to the particle's rest-frame. In it was shown 
how to extend the standard variational formulation of field theory to the case of non-inertial 
frames (for the definition of the co- moving frame see also 0). In the present Section we 
show how to extend the standard Hamiltonian formulation of a classical field theory. 

Consider any field theory based on a first-order relativistically-invariant Lagrangian den- 
sity 

L = L{i;,di;), (3.1) 

where ip is a (possibly multi-index) field variable, which we do not need to specify at the 
moment. As an example, ip could denote a scalar, a spinor or a tensor field. 

We will describe the above field theory with respect to accelerated reference frames, 
related with observers moving along arbitrary space-time trajectories. Let ( be such a 
(time-like) trajectory, describing the motion of our observer. Let = q^{t), k = 1,2,3; be 
the description of ( with respect to a laboratory reference frame, i. e. to a system (y'^), A = 
0,1,2,3; of Lorentzian space-time coordinates. We will construct an accelerated reference 
frame, co-moving with (. For this purpose let us consider at each point (t, q(t)) G ( the 
3-dimensional hyperplane Ht orthogonal to (, i.e. orthogonal to the four-velocity vector 
U{t) = (u'^it)): 

{un = {u',u''):=j{l,v') , (3.2) 

where := q^ and the relativistic factor 7 := 1/Vl — v^. We shall call the "rest frame 
surface". Choose on any system (x*) of cartesian coordinates, such that the particle is 
located at its origin (i. e. at the point x' = 0). 
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Let us consider space-time as a disjoint sum of rest frame surfaces S^, each of them 
corresponding to a specific value of the coordinate := t and parameterized by the coor- 
dinates (x*). This way we obtain a system (x") = {x^,x^) of "co-moving" coordinates in a 
neighbourhood of (. Unfortunately, it is not always a global system because different S's 
may intersect. Nevertheless, we will use it globally to describe the evolution of the field ip 
from one to another. For a hyperbolic field theory, initial data on one St imply the entire 
field evolution. We are allowed, therefore, to describe this evolution as a 1-parameter family 
of field initial data over subsequent S's. 

Formally, we will proceed as follows. We consider an abstract space-time M := T x S 
defined as the product of an abstract time axis T = with an abstract, three dimensional 
Euclidean space S = R^. Given a world-line (, we will need an identification of points of M 
with points of physical space-time M. Such an identification is not unique because on each 
Sj we have still the freedom of an 0(3)-rotation. 

Suppose, therefore, that an identification F has been chosen, which is local with respect 
to the observer's trajectory. By locality we mean that, given the position and the velocity 
of the observer at the time t, the isometry 

F(q(t),v(t)) : S 1-^ St (3.3) 

is already defined, which maps G S into the particle position {t, q(t)) G S^. 

As an example of such an isometry which is local with respect to the trajectory we 
could take the one obtained as follows. Choose the unique boost transformation relating the 
laboratory time axis d/dy^ with the observer's proper time axis U. Next, define the position 
of the d/dx'^ - axis on St by transforming the corresponding d/dy^ - axis of the laboratory 
frame by the same boost. It is easy to check, that the resulting formula for F reads: 

y\t,x') := t + -i{t)x'vi{t) , 

y\t,x') := q\t) + {51 + ipiw^yvi) x' . (3.4) 
Here, the following function of a real variable has been used: 

The function is well defined and regular (even analytic) for = r < 1. The operator 

^-^a'^^ := 5f + i^v\ (3.6) 

acting on rest-frame variables x' comes from the boost transformation (for simplicity we skip 
the argument of the function ^p). 

Suppose, therefore, that for a given trajectory C, a local isometry (|3.3|) has been chosen, 
which defines : M i— > M. This mapping is usually not invertible: different points of M 
may correspond to the same point of space-time M because different St's may intersect. It 
enables us, however, to define the metric tensor on M as the pull-back F^g of the Minkowski 
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metric. The components ga/B of the above metric are defined by the derivatives of F^, i. e. they 
depend upon the first and the second derivatives of the position q(t) of our observer. 

Because [x'^) are cartesian coordinates on E, the space-space components of g are triviah 
9ij ~ ^ij- '^^^ only non-trivial components of g are, therefore, the lapse function and the 
(purely rotational) shift vector: 

N = -yi=^ = + aix') , 

Nm = gOm = J~'^emklUj''x'- , (3.7) 

where a* is the observer's acceleration vector in the co-moving frame. The quantity uj"^ is 
a rotation, which depends upon the coordination of isometries ( |3.3| ) between different S^'s. 
Because u"^ depends locally upon the trajectory, it may also be calculated in terms of the 
velocity and the acceleration of the observer, once the identification ( |3.3| ) has been chosen. 
In the case of example (|3.4|) , it is easy to check that 

d = 8i\v' , (3.8) 



f^m = l^mlV^ , (3.9) 

where 

<^m« := l^v'^eum , (3.10) 

and i)^ is the observer's acceleration in the laboratory frame. 

The metric F^g is degenerate at singular points of the identification map (i.e. where the 
identification is locally non- invert ible because adjacent S's intersect, i. e. where iV = 0), but 
this degeneration does not produce any difficulties in what follows. 

3.2 Lagrango-Hamiltonian 



Once we know the metric (|3.7|) on M, we may rewrite the invariant Lagrangian density L of 
the field theory under consideration, just as in any other curvilinear system of coordinates. 
The Lagrangian obtained this way depends on the field '0, its first derivatives, but also on 
the observer's position, velocity and acceleration. Variation with respect to ip produces field 
equations in the co-moving system (x"). Due to the relativistic invariance of the theory, 
variation of the Lagrangian with respect to the observer's position q should not produce 
independent equations but only conservation laws, implied already by the field equations. 

For our purposes we will keep, however, at the same footing the field degrees of freedom 
if) and (at the moment, physically irrelevant) observer's degrees of freedom q^. We are going 
to perform the complete "Hamiltonization" of this theory, i.e. to pass to the Hamiltonian 
description both in field and observer's variables. 
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Let us first perform a partial Legendre transformation, and pass to the Hamiltonian 
description of the field degrees of freedom, keeping the Lagrangian description of the "me- 
chanical" degrees of freedom. For this purpose we define 

LH:=L-n^, (3.11) 

where 11 is the momentum canonically conjugate to ip: 

n:=^. (3.12) 

The function Lh plays the role of a Hamiltonian (with negative sign) for the fields and a 
Lagrangian for the observer's position q. It is an analog of the Routhian function in analytical 
mechanics (cf. 0). The Lagrango-Hamiltonian generates the Hamiltonian field evolution 
with respect to the accelerated frame, when the "mechanical degrees of freedom" are fixed. 



Due to (|3.7| ), this evolution is a superposition of the following three transformations: 

• time-translation in the direction of the local time-axis of the observer, 

• boost in the direction of the acceleration of the observer, 

• purely spatial 0(3)-rotation uj^. 

It is, therefore, obvious that the numerical value of the generator of such an evolution is 
equal to 

Lh = {n + a'TZk - u^Sm) , (3.13) 

where 7i is the rest-frame field energy, TZk is the rest-frame static moment and Sm is the 
rest-frame angular momentum, all of them calculated on S. The factor 7"^ in front of the 
generator is necessary, because the time t = x^, which we used to parameterize the observer's 
trajectory, is not the proper time along ( but the laboratory time. 

3.3 Legendre transformation and relativistic invariance 

Now, we perform the Legendre transformation also with respect to the observer's degrees of 
freedom, and find this way the complete Hamiltonian of the entire (observer + field) system. 
Let us observe that Lh is a 2-nd order Lagrangian in the observer's variable: 



LH(q,v,V, fields) , (3.14) 



(for the Hamiltonian description of a theory arising from a 2-nd order Lagrangian see Ap- 
pendix A general discussion of a Hamiltonian formalism arising from higher order La- 
grangians may be found in |16[ - ||18|). Let pk and vTfc denote the momenta canonically 



conjugated to q and v respectively. The phase space of the entire system, parameterized 
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by {q'^yv'^yPkyTik) in the observer's sector and by {ip,H) in the field sector, is endowed with 
the canonical symplectic 2-form: 

Q := dpk A dq'' + dTik A dv'' + Q^'^^'^ , (3.15) 

which generates the Poisson bracket for any two observables and Q: 

^ '^^ ■ [dq^dpk dpkdq^) \dv^d-nk ^77^9^;^ ' ' 

Due to the relativistic invariance of the theory the following observables: Ti, TZk, Sm and 
the rest-frame field momentum Vk generate with respect to (|3.16|) the Poincare algebra: 





0, 




0, 




-Vk , 




0, 


{Vk^niY'^'^ = 


—Qki^ , 


[Vk^SiY'^"' = 


(-klrnP ) 


{TZk^niy'''" = 


~^klmS 


{TZk^Siy'^'' = 


^klrn^ 1 


{Sk,Siy''''' = 


am 
tklm'J 



(3.17) 



To perform the Legendre transformation one has to calculate q^ and v'' in terms of q^, ^ pk 
and TTfc from formulae (cf. Appendix]^): 

= ITT ' '^k = ITT ■ ^-^^ 
Since Lh is linear in v the Legendre transformation is singular and it gives rise to the 



Hamiltonian theory with constraints (see ||T^- ||14|). The primary constraints follow from 
iH): 

0^1^ := VTfc - ^ = TTfe + (a^,7^^ - u^^kS"^) ^ . (3.19) 

By the weak equality symbol "~" we emphasize that the quantity 0*'^^ is numerically re- 
stricted to be zero but does not identically vanish throughout the phase space. This means 
in particular, that it has nonzero Poisson bracket with the canonical variables. This way the 
complete Hamiltonian reads 

H = pkV^ + vTfci;^ -Lh= Pkv' + i-'n + i^V^fc • (3.20) 



The observer's acceleration v plays in ( |3.2CI| ) the role of Lagrangian multiplier. To reduce the 
theory with respect to the constraints (|3.19| ) and to find the "true" degrees of freedom one 
has to find the complete hierarchy of constraints via so called Dirac-Bergmann procedure 
(cf. 0-0). 
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PROPOSITION 1 The secondary constraints cj) i! read: 

cf>^'^ := {^^Z H} = -Pk + IVkH + ^ . (3.21) 

Let us find the constraint algebra, i.e. "commutation relations" between constraints 0^^''. 
PROPOSITION 2 

{0(J,0?} = O, (3.22) 

for a,b = 1,2 and k,l = 1, 2, 3. 

It means that (p^]^ and 0*-^^ are so-called first-class constraints (see [|l^ - [1^). One has to 
continue this algorithm and look for the constraints which are the "conservation laws" for 
the 0^J, i.e. 



k 



(=')-{0(J,H}^O. (3.23) 

(3) 

However, the tetriary constraints ^ are satisfied identically due to the Poincare relations 



This way we show that the hierarchy of constraints ends on the level of secondary con- 
straints and we obtain the Hamiltonian theory with 6 first-class constraints 0*"^^ Let V 
denotes the constraint subspace of the phase space V, i.e. 

r-={xeP I 0^J (a;) = for a = 1, 2; k = 1, 2, 3} , 

and let e : P P be an embedding. Then, due to ( p.22| ), the pull-back e*fl of the symplectic 
form f2 is degenerate and to each functional 0*'^^ corresponds a "gauge direction" X^'^ , i.e. 

J e*Q = , 

such that 

e,x(^)jfi = 50(J , 

i.e. X^l^ is a Hamiltonian vector field corresponding to the functionial 0*'^^ {X^l^ is tangent 
to V due to (13:2^)). 

Therefore, our theory is a gauge theory with 6 gauge parameters: q and v. Canonically 
conjugated momenta p and tt are subjected to constraints 0^^''. Thus, reducing the theory 
with respect to these constraints (i.e. passing to gauge equivalence classes) we end up with 
the "true" degrees of freedom, namely those describing the field. Fixing the trajectory plays 
the role of "gauge fixing" and the "evolution equations" of the observer are automatically 
satisfied if the field equations are satisfied. 

This result is an obvious consequence of the relativistic invariance of the theory. Ob- 
server's parameters are not true degrees of freedom and can be easily "gauged away". 
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4 Renormalized electro dynamical Routhian 



We would like to apply the formalism presented in the previous Section for the Hamiltonian 
description of electrodynamics of a point particle. It turns out that the use of the particle's 
rest-frame simplifies considerably the formulation of the theory. Hence we will use not an 
arbitrary observer, but the one following exactly the particle's trajectory. As a consequence, 
the trajectory will no longer be a gauge parameter, but will have an independent, dynamical 
meaning, as a new degree of freedom of the theory. 

Let us observe that starting from the Lagrangian descripton of electrodynamics given by 
the standard Maxwell Lagrangian, the Legendre transformation ( |3.11| ) does not lead to the 
correct local expression for the field energy (we obtain the "canonical Hamiltonian" which 
differs from the field energy by a complete divergence). However, in the "field sector" we 
may already start with the electrodynamical Routhian (Lagrango-Hamiltonian) (|3.13|) . To 
obtain this generator we take as 7^, TZk and Sm the conventional energy, static moment and 
angular momentum of the electromagnetic field. These quantities are defined as appropriate 
integrals of the components of the symmetric energy-momentum tensor: 



In the Section ^ we construct the Hamiltonian structure of Maxwell electrodynamics which 
is perfectly suited for the formalism introduced in the previous Section. In particular, the 
"new Hamiltonian" is positive definite and equals to the field energy, i.e. is defined as an 
integral of T'^. 

When one adds a point particle to the electromagnetic field, then, obviously, the total field 
energy is not well defined due to the singularity of the particle's Coulomb field. Therefore, 
one has to perform renormalization. Decomposing the electric induction field on the rest- 
frame surface into the sum 

D = Do + D 

of the Coulomb field Dq = and the remaining part D, we obtain the following formulae 
for the renormalized rest-frame quantities (see and ^ for details): 

n = m + ^J^CD' + B^)d'x, (4.1) 

Vi = / (D X B)i £x + [{Box B)i Sx , (4.2) 

TZk = \ I Xfc(D' + B2)rf3x+ / Xkmod'x , (4.3) 

Sm = [ emkix^BxByd'x . (4.4) 

Finally, the renormalized Lagrango-Hamiltonian Lh is defined as follows: replace in ( |3.13|) 
the exact values, i. e. calculating for the complete field theory, of the quantities Ti, TZ^ and 
Sm by the above renormalized quantities (f4.1|) , (4.3) and (|4.4| ), containing only the external 
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Maxwell field. It was shown in that the variational principle applied to the renormalized 
Lh gives the Euler-Lagrange equations which are equivalent to the fundamental equation 



5 The new gauge- invariant Hamiltonian structure of 
Maxwell electrodynamics 

In field theory, contrary to the classical mechanics, there is no unique way to represent the 



field evolution as an infinite-dimensional Hamiltonian system (see [10| and Each such 

representation is based on a specific choice of controlling the boundary value of the field, and 
corresponds to a specific choice of the Hamiltonian. This non-uniqueness is implied by the 
non-uniqueness of the evolution of the portion of the field, contained in a finite laboratory 
V . Indeed, the evolution is not unique because external devices may influence the field 
through the open windows of our laboratory. To choose the Hamiltonian uniquely, we have 
to insulate the laboratory or, at least, to specify the influence of the external world on it. 
One may easily imagine an unsuccessful insulation, which does not prevent the external field 
from penetrating the laboratory. From our point of view, an insulation is sufficient if it keeps 
under control a complete set of field data on the boundary dV in such a way, that the field 
evolution becomes mathematically unique. 

For relatively simple theories (e.g. scalar field theory) the Dirichlet problem may be 
treated as a privileged one among all possible mixed (initial value + boundary value) prob- 
lems which are well posed. It turns out (see |^) that the Dirichlet problem leads to the 
positive definite Hamiltonian. This means that there is a natural way to insulate the labo- 
ratory V adiabatically from the external world. But already in electrodynamics (and even 
more in General Relativity) any attempt to define the field Hamiltonian leads immediately 
to the question: how do we really define our Hamiltonian system? 

5.1 Canonical approach 

One usually starts with a canonical dynamical formula for electrodynamics (see 0, |ll|): 

- ifc(5J^'^° = -5/7^?"+/ J^'^^M, , (5.1) 

V JdV 

where Hy"' is a "canonical Hamiltonian" related via the Legendre transformation to an 
electro dynamical Lagrangian density C : 

Tcan I I -7-fcO 



We use standard notation: J-'^'^ denotes the electromagnetic induction tenor-density. The 
volume V belongs to the hyperplane and consists of the exterior of the sphere 5'(ro) (by 
_L we denote the component orthogonal to the boundary). We stress that the formula ( |5.1| ) 
is true also in the case of nonlinear electrodynamics. For the linear Maxwell theory C is 
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the standard Maxwell Lagrangian L Maxwell = —\^^'^fij,u, where /^,^ := d^A^ — d^A^ is the 
electromagnetic field tensor. In this case the electromagnetic induction tensor-density is 
given by J^^'' := ^/^g^'^g"'^ fap- 

To describe the boundary term it is convenient to use spherical coordinates (^'^), a = 
1,2,3; adapted to dV. We choose = r as the radial coordinate and (C,^), A = 1,2; as 
angular coordinates: = 9, = ip. The Euclidean metric gab is diagonal: 

6-33 = 1 , gu = r'^ , 6-22 = sin^ 6 , 
and the volume element A = (det gabY^"^ is equal to sin B. With this notation we have: 

/ J^*^°5Afc - AkSJ^''^ = + / ^°=^Mo + J^^'^SAb . (5.3) 

Jv JdV 

The formula (|5.3| ) is analogous to the Hamiltonian formula in classical mechanics 

Pkdq'' - q^dpk = -dH{q,p) . 



But there is also a boundary term in ( p.3| ), typical for field theory. Killing this term by an 
appropriate choice of boundary conditions is necessary for transforming the field theory into 



an (infinite dimensional) dynamical system (see 0, [10 1). Thus, boundary conditions for 
y4o|c?\^ and A^I^T^ make the electrodynamics equivalent to the infinite-dimensional Hamil- 
tonian system. From the mathematical point of view this is the missing part of the definition 
of the functional space. The Hamiltonian structure of Maxwell electrodynamics described 
above is mathematically well defined, i.e. a mixed Cauchy problem (Cauchy data given on 
and Dirichlet data given on dV x R) has an unique solution (modulo gauge tranformations 
which reduce to the identity on dV x R). 

5.2 The new approach 

It turns out that there is another Hamiltonian structure which is also mathematically well 



defined (see [|rT], [|T^). Moreover, it possesses very nice properties from the physical point 
of view: 

1. the new Hamiltonian corresponds to the field energy obtained from the symmetric 
energy-momentum tensor, i.e |(-D^ + -B^) in the laboratory frame, which, contrary to 
the canonical Hamiltonian, is positive definite, 

2. this structure is perfectly suited for the reduction of the theory with respect to the 
Hamiltonian constraint D^,k = , i.e. Gauss law (see [n|). 



Let us start from the canonical relation (|5.3|) 

= + / j^^^As + j^^^As - A^5J^^^ - is^J^^o . (5.4) 

Jv 



L 
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On each sphere S{r) = {r = const} the 2-dimensional covector field Ab sphts into a sum of 
the "longitudinal" and the "transversal" part: 



Ab = u,B+eB'^v,c , (5.5) 

where the coma denotes partial differentiation and e"^^ is a sqew-symmetric tensor, such that 
\e^^ is equal to the Levi-Civita tensor-density (i. e. Ae^^ = — Ae^^ = 1). The functions u 
and V are uniquely given by the field up to additive constants on each sphere separately. 
Inserting this decomposition into ( |5.4]) and integrating by parts we obtain: 



Jv 

= + / J^°3(5Ao - J'^'^B 5u - J^^^c ej^6v . (5.6) 

JdV 

Using identities DbT^^ + c^sJF^^ = and BbT^^ + d^T^^ = 0, implied by the field equations 
d^T^^ = 0, and integrating again by parts we finally obtain: 

/ [J^''6{A, - u,3 ) - (is - ^,3 )5^'°] + [-^°^ll^6^c5t; - vdT'^^^'^eBc] 

= + / ^°3(5(Ao -u)+ ^^^ll^eBC^^ • (5.7) 

JdV 

Here, by "||"we denote the 2-dimensional covariant derivative on each sphere S{r). The 
quantities {Aq — u,o ) and (^3 — u,^ ) are "almost" gauge invariant: only their monopole 
part (mean-value) on each sphere may be affected if we change the additive constant in the 
definition of u (the choice of an additive constant in the definition of v is irrelevant, because 
it is always multiplied by quantities which vanish when integrated over a sphere). The sum 



of the volume and surface integrals in (|5.7]) is however gauge invariant. Now, 

= (curlA)^ = e^^(AB,3 -^3,5 ) (5.8) 

and using ( |5.5| ) we have 

A{As-us)=r'B^\\''eAB, (5.9) 

where A denotes the 2-dimensional Laplace-Beltrami operator on S{r) multiplied by (the 
operator A does not depend on r and is equal to the Laplace-Beltrami operator on the 
unit sphere 5'(1)). The operator A is invertible on the space of monopole-free functions 
(functions with vanishing mean value on each S{r)). This functional space will play an 
important role in further considerations and all the dynamical field quantities of the theory 
will belong to this space. To fix both terms in ( ^.7| ) uniquely we choose u in such a way that 
the mean value of (^3 — ^^3) vanishes on each sphere. Hence, with the above choice of the 
additive constants the quantity ^3 — u 3 becomes gauge invariant: 

A,-u, = r'A-\B^\\^eAB). (5.10) 
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Let us observe that the function v is also gauge invariant (up to an additive constant, 
which does not play any role and may also be chosen in such a way that its mean value 
vanishes on each sphere). Indeed, we have: 

= (curl A)=^ = A^^^^esA = -r'^Aw . (5.11) 

Due to the Maxwell equation div B = 0, the function is monopole-free and the Laplasian 
A may again be inverted: 

v = -r^A-^B^. (5.12) 



The formula (|5.7|) could be also obtained directly from (|5.4| ) by imposing the following gauge 
conditions: 

= 0, (5.13) 
AAs = 0. (5.14) 

Sir) 

The above condition does not fix the gauge uniquely: we still may add to A^ the gradient 
of a function of time / = f{t). This residual gauge changes only the monopole part of Ao, 
but both the volume and the surface integrals in (|5.3| ) remain invariant with respect to such 
a transformation. 

Assuming the above gauge, we have u = and As = e^^v,c- To simplify the notation 
we will, therefore, replace our invariants {Aq — u,o ) and (A3 — ^,3 ) by the values of Aq and 
A3, calculated in this particular gauge. 

Let us observe that formula (|5.7|) represents dynamical system with infinitely many de- 
grees of freedom described by four functions: A3, J-'^^^^'"eBC and v (contrary to ( |5.3| ) 
described by six functions A^ and D^). Two of them will play the role of field configurations 
and the remaining two will be the conjugate momenta. Let us consider a boundary term 
in ( |5.7| ). Killing this term by an appropriate choice of boundary conditions is necessary for 
transforming the field theory into an (infinite dimensional) dynamical system (see |^, [|T0[]). 
/^From this point of view, the quantity v (or, equivalently B^) is a good candidate for the 
field configuration, since controlling it at the boundary will kill the term 6v in the boundary 
integral. On the contrary, 6Aq can not be killed by any simple boundary condition imposed 
on A3. We conclude, that it is rather J-'^^ = XD^ than A3, which has to be chosen as another 
field configuration. 

Hence, we perform the Legendre transformation in formula (|5.7|) on the boundary dV: 
J^^^SAo = 5(J^°Mo) - Ao6T^\ This way, using (|ig) and (|1|), we obtain from (|3) the 
following result: 

/ Ar4 - Z)35(A-i(5^ll Vb)) + A-\B^^^^eAB)SD' 
Jv 

+ BH{A-\D''\\''eAB)) - A-i(I)^ll%B)5i?' 
= sfm^''- [ ^03^0)+/ -XAo5D^ -r^A-\j^^^\\''eAB)SB'\ (5.15) 

V JdV / JdV 
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where we used the fact that the operator is self-adjoint on the functional space of 
monopole-free functions on a sphere. 

We see that {D^, B^) play the role of field configurations, whereas the remaining functions 
[B'^^^^tABi -D^"^e^B) describe the conjugate momenta. Controlling the configurations at the 
boundary we kill the surface integral over dV and obtain this way an infinite dimensional 
Hamiltonian system describing the field evolution. There is, however, a problem with such a 
control, because the electric induction cannot be controlled freely on the boundary. The 
reason is that the total electric flux through both components of dV (i. e. through S'(ro) and 
through the sphere at infinity) must be the same: 

5(ro) Js{rac) 

where e is the electric charge contained in S'(ro). Hence, we have to separate the monopole- 
free ( "radiative" ) part of (which can be independently controlled on both ends of V) from 
the information about the electric charge. For this purpose we split the electric induction 
into 

D' = -^^+^\ (5.17) 



where D is again a monopole-free function. It follows from (|5.16|) that the monopole part 



of (equal to e/47rr^) is nondynamical and drops out from the volume integral of 
because it is multiplied by a monopole-free function B^^^^tAB- The remaining part (which 
does not carry any information about the charge e), together with B^, can be taken as the 
true, unconstrained degrees of freedom of the electromagnetic field. 
In the same way we split the scalar potential Aq. 

Ao = + Ao, (5.18) 

where 0(r) is the mean value of Aq on the sphere S{r) (monopole part) and Aq is a monopole- 
free function ("radiative" part of Aq). Now, the boundary term Aq5D^ in ( |5.4|) reads 

XAo6D^ = ^f Xr-^(j)6e+ f XAo6d\ (5.19) 
av 47r Jdv Jav 

Finally, we perform the Legendre transformation between and the monopole part of at 
infinity. Hence, we control the total charge contained in S'(ro) and the monopole function cj) 
at infinity. Since the latter does not contain any physical information and is used only to fix 
the residual gauge, we may use the simplest possible choice: 0(oo) = 0. This way we have 
proved the following 



THEOREM 3 The quantities 



^1 = rS^ 

*2 = rlf , 

XI = -rA-i(Z}-^IIVB) 

X2 = r/\-\B^\\^eAB) 
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together with the value e of the electric charge contained in S'(ro) contain the entire (gauge 
invariant) information about the electromagnetic field. 

Quantities play the role of field configurations and xa are conjugated momenta. Finally, 
formula ( |5.15| ) reads: 

/ X{xaS^^ - ^^Sxa) = -SHv + I \Xa^^^ + 6^0(00) + 0(ro)(5e (5.20) 
Jv Jav 

where the new Hamiltonian Hy equals 

Hv = - / ^°'Ao , (5.21) 

JdV 

and the "boundary momenta" are given by: 

X2 = -r-'Ao. 

They describe the response of the system to the control of the boundary values of configura- 
tions The above Hamiltonian corresponds to the symmetric energy-momentum tensor 
of the Maxwell field (cf. IT^, 0]), i.e. Hy equals numerically to the amount of the elec- 
tromagnetic energy contained in the volume V. Obviously, the limit lim^o^o Hy is not well 
defined due to the Coulomb field singularity. However, renormalizing it one gets exactly (up 
to a sign) the renormalized Lagrango-Hamiltonian Lh- 
The space of the electromagnetic field contained in V 

pfieid ^ i^A^ ^^,Y with boundary conditions ^"^{dV} 

is endowed with the canonical symplectic 2-form 

^f ieid ._ ^^^^ ^ ^^A _ ^5 22) 

It may be easily obtained (see |Ty]) by the reduction of the standard presymplectic form 
Jy A with respect to the Hamiltonian constraint d^J^^^ = 0. The form flf.^'^'''^ defines 
in the space of physical observables, i.e. functionals over P/j'^''^, the canonical Poisson bracket 
structure: 

{T Cl/-"^ ■= f A ( - _^JL\ (K 9.1 

^ ■ U^^(x)5xa(x) 5^^i^) 6xAi^) ) ■ ^ ^ 



6 Reduction 

In this Section we finally find the space of Hamiltonian variables and the total Hamilto- 
nian for electrodynamics of a point particle. We start with the renormalized Lagrango- 
Hamiltonian. It generates the Hamiltonian dynamics in the "field sector" and Lagrangian 
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dynamics in the "particle's sector". Our aim is to perform complete "Hamiltonization" of 
the theory, i.e. to perform the Legendre transformation in the particle's variables. 

The original phase space is a direct sum of a particle's space pp"''**^^*^ and a phase space 
of the Maxwell field V^'^^'^: 



yjj yjjparticle field 



5 



where 



V^'"^"^ := lim P/;'^''^ . 



The particle's phase space pp'*''**'^^^ is parameterized by q'^, and conjugated momenta pfc, 
TTfc and T'f'^^'-'^ by and conjugated momenta xa with the boundary condition (now 
V is an exterior of the sphere S{ro) with rg ^ 0). The total phase space V is endowed with 
the canonical symplectic form Q given by 

n := dpk A dq'' + diTk A dv'' + (6.1) 

which generates the Poisson bracket in the space of functionals over V: 

{TM := f^^-^^Uf— (6.2) 

where 

n^'^^"^ ■= lim nli^^"^ , (6.3) 

ro->0 

{-p^gy^eid linil^,^}/-"^ , (6.4) 

and { , is given by ( |5.23| ). In the case of complete field theory the rest frame functionals: 
H, Vk, Tlk and Sm form the Poincare algebra with respect to { , It is no longer true 

for the renormalized functionals. This does not mean that renormalized electrodynamics is 
incompatible with the special theory of relativity. Notice, that already in the context of 
inhomogeneous Maxwell theory with given point-like sources, field functionals do not form 
the Poincare algebra but the theory is obviously relativistically invariant. We have shown 
in Section |] that the Poincare algebra structure is equivalent to the fact that the Legendre 
transformation ( |3.11D leads to the Hamiltonian theory with first class constraints. A theory 
which is subjected to the first class constraints is a gauge-type theory. Now, we expect 
that the particle's trajectory is no longer a gauge parameter but plays a dynamical role. 
Therefore, the "breaking" of the Poincare algebra structure is the necessary condition for 
the nontrivial particle's dynamics. 

To effectively reduced V we shall proceed as follows: we consider particle's trajectory C, 
as a limiting case of a tiny world-tube of radius tq. Therefore, calculating Poisson brackets 



according to ( |6.2| ) we shall keep tq > and then finally go to the limit tq ^ 0. 
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One can easily shows that only three among nine Poincare relations (|3.18|) are changed, 
namely: 

{Vk,Vi}t^'' = eumB'^iro), (6.5) 



QirrQ J 

where B"^{ro) denotes the mean value of the m-th component of B on the sphere 5'^(ro). 
Functions and Pk are components of the dipole parts of the functions a and /?, given by 
decompositions: 

DP(a) = aky , DP(/3) = A— . 

Functions a and f3 are and r° terms respectively in the decomposition of the radial 
component in the vicinity of the particle (see ( p.6| )). 

Let us perform the Legendre transformation in the particle's sector in the same way as in 
the Section ^. The renormalized Lagrango-Hamiltonian gives rise to the Hamiltonian theory 
with constraints. The primary constraints (j)^]} are given obviously by (|3.19|) . It is easy to 
show that the secondary constraints 0^ 

0^? := , (6.6) 

where the complete Hamiltonian H is given by ( |3.20| ), have the same form as in ( p.21[ ). The 
next step in the Dirac-Bergmann procedure is to calculate the Poisson brackets between 
constraints. 

PROPOSITION 3 

{0^l\0^P}ro = 0, 

0^?}ro = j[m- J [gu + l^VkV^ , (6.7) 



0, 



{0^iU^?}ro = 27eV^^](£jr + A)+ee,,„5"^(ro) 



■where 1A\^kBi\ := A^Bi - AiBk- 



i^i^ and 

theory, however, the constraint algebra is completely different. 



We conclude that the constraints 0^ and (f) p! are the same as in the case of fundamental 



DEFINITION 2 A functional T over V is said to be first-class if its Poisson bracket with 
every constraint vanishes weakly. 
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Observe that due to ( |6.7|) there is no functional among 0*'^-' which is first-class. The absence 
of the first-class constraints means that there are no gauge parameters in the theory. The 
constraints which are not first-class are called second-class. We conclude that the renormal- 
ized Lagrango-Hamiltonian gives rise to the Hamiltonian theory with second-class constraints 

Let V denotes the constrained submanifold of V, i.e. 

P:={xGP|0(J(x)=O}, 

and let e : P ^ P be an embedding. Due to the fact that 0^^'' are second-class constraints 
the pull-back Q = e*Q of the symplectic form f2 is a non-degenerate 2-form on V (this is 
equivalent to the fact that there is no gauge freedom at all). One may think that {V, Q) 
is an adequate phase space for electrodynamics of a point particle. But it is not the whole 
story. 

PROPOSITION 4 In the finite- dimensional case a symplectic form defines an isomor- 
phism between vectors and covectors (1-forms), i.e. if {P,uj) is a symplectic manifold, then 
ijj induces a continuous linear map for any point p G P: : TpP T*P defined by 

ujliX)-Y ■.= Up{X,Y) , 

for any X,Y G TpP. If dimP < oo, then cUp is an isomorphism. 

This theorem, however, is no longer true in the infinite-dimensional case. 

PROPOSITION 5 If dim P = oo, then in generalcUp is only injective, i.e. iftUp{X,Y) = 
for all Y e TpP, then X = 0. 



See ||I5[ for the proof. 

DEFINITION 3 If is only injective, then u is called a weak symplectic form. If, more- 
over, iUp is onto, then cu is called a strong symplectic. 

Let (P, u) be a weak symplectic manifold and let X be a vector field on P defined on a dense 
subset D of P. 

DEFINITION 4 We call X a Hamiltonian vector field if there exists a functional JF : D — > 
R such that 

X_\lo = -5T 

is satisfied on D. 
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Let us observe, that for a weak symplectic u, there need not exist a vector field Xyr cor- 
responding to every given functional J-" on D. Therefore, from the physical point of view, 
weak structure is "too weak" . 

Let us come back to a manifold V with a nondegerate 2-form Q. Is {V, fl) weak or strong? 
It turns out that Q is only a weak symplectic form on V and the reason is very physical. 
The Hamiltonian vector field corresponding to the Hamiltonian H is not well defined on 
V. To simplify our considerations let us parametrize P by a suitable coordinate system. 
The simplest parametrization is the following: field sector is still parameterized by \E'^ and 
Xa and particle's sector by q and v. Momenta pk and vr^ are completely determined by 
constraints (|319D and ( ^ ): 

7rfc(q,v, fields) = —f~\a},ni-u;^kS"'), (6.8) 
Pfc(q,v, fields) = ^Vkn + ^-^8L\Vi . (6.9) 

In fact TTfc and pk do not depend on q. Notice, that momentum p^ canonically conjugated 
to particle's position equals to the total momentum of the composed system (particle + 
field) in the laboratory frame. Now, using ( p.20|) and (|6.9| ) the complete Hamiltonian H on 
V is given by 

H(q, v; fields) = 1 [H + v'^Vk) , (6.10) 

and equals to the total energy of the composed system in the laboratory frame. Obviously, 
H is well defined at any point of V. However, the Hamiltonian vector field Xu is defined 
only on a subset V* of V. 

THEOREM 4 The Hamiltonian vector field Xn is well defined if and only if the 'funda- 



mental equation" (WW is satisfied, i.e. V* is defined by the following condition: 



DP(47rm^2 + e^^^g )(0) = . 

The proof of Theorem ^ is given in the Appendix The main result of this Section consists 
in the following 

THEOREM 5 {V*,Q*) is the strong symplectic manifold. 

The proof of this Theorem is given in the next Section where we construct the reduced Poisson 
bracket on V* . It turns out that the reduced bracket of any two well defined functionals over 
V* is well defined throughout the reduced phase space. 

Therefore, we finally take the space {V*,Q*), where Q* is a reduction of Q to V*, as 
a phase space for the composed system. Due to Theorem | the Hamiltonian structure 
{V*, Q*, H) for electrodynamics of a point particle is well defined: each state in V* uniquely 
determines the entire history of the system. Observe, that due to ( |6.8| ) and ( |6.9| ) the reduced 
symplectic form possesses highly nontrivial form: 

n* = ^dv' A dq' + / X^J^!^ Adq'^+ f X^Sxa A dq' 



23 



+ / A5xa A5^^ . (6.11) 

It was shown in |0 that the "fundamental equation" is equivalent to the total momentum 
conservation. In the Lagrangian formulation of the theory [0 this equation is nothing else 
that the Euler-Lagrange equation of the variational problem. Now, in the Hamiltonian 
formulation, the "fundamental equation" is already present in the definition of the phase 
space of the system and the name "fundamental" is fully justified. 



7 The Poisson bracket 

In this Section we find the Poisson bracket structure for electrodynamics of a point particle, 
i.e. we find the reduced Poisson bracket { , }* for the functionals over V*. It is defined in a 
obvious way via the symplectic form Q* on V*: 

DEFINITION 5 For any two functionals T and Q over V* 

To find the explicite form of { , }* let us apply the Dirac method [|12| (that is way some 
authors, e.g. [|l^ - [Q, call it the Dirac bracket). As in the previous Section we find it 
firstly for rg > and then go to the limit tq ^ 0. 

Using constraint algebra ( |6.7|) -( |6l8D let us define a following matrix (so-called Dirac ma- 
trix): 

Notice, that on V* , i.e. on a dense subset of V where the "fundamental equation" is satisfied, 
the constraint algebra reduces to 

Xijiro) = 1 - {gij + I'^ViVj^ , (7.1) 

y,,(ro) = 2-^V[,(3,] L- +ee,,„S'"(ro) . (7.2) 

m \ OTrro / 

For a Hamiltonian theory with second-class constraints the Dirac matrix is non-singular and 
the inverse to Cjj(ro) reads: 

^ - [ A(ro) j ' 
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where 



A^^iro) := {X-\ro)r = ^ , \ (7.3) 



m — t; 

on ro 



B'^iro) := {X-\ro)Y{ro)X{ror'y^ 



mm—-^— (m — TT-— )^ 

Finally, the "Dirac bracket" is defined as follows 

{^,^r:= hm{^,^};, (7.5) 

where 

e;}; := gU + ^''(^O) ({-^, 0?}ro{0^ ^In. - 0?}ro{0^;\ ^}ro^ 

- 5^^(ro){^,0(l)U{0^;\^}ro. (7.6) 

We stress that the reduced Poisson bracket { , }* is well defined for any two well defined 
functional T and Q on the reduced phase space V* . Moreover, particle's and field degrees 
of freedom are kept at the same footing. To our knowledge it is the first consistent Poisson 
bracket structure for the theory of interacting particles and fields. 

The complete set of "commutation relations" defined by { , }* is given in the Appendix 0. 
Using these relations one easily prove the following 

THEOREM 6 Hamilton equations: 

|vl;^(x) = {vI/^(x),Hr, 

|xa(x) = {xA(x),Hr, 

reconstruct Maxwell equations in the co-moving frame (cf. j^). Moreover, 

q' = {q',lir = v' , 

and 

t'' = {t^',Hr = -(a-i)'=,A , 
m 



where the inverse of the operator a defined in ( \3.(^ reads 

gives the particle's "equation of motion": 

mak = e(3k . 
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8 The Hamiltonian 



The formula ( |6.1(j| ) gives the quasi-local Hamiltonian for the composed "particle + field" 
system. It is expressed in therms of q, v and field variables. However, these variables 
are highly noncanonical with respect to { , }*. Let us observe that momentum p and the 
particle's position q are still conjugated to each other with respect to the reduced bracket 
(|7.5| ). Indeed, one easily shows that 

Therefore, we can use p instead of v to parametrize V* . However, to express the complete 
Hamiltonian H in terms of p one has to express v in terms of q, p and fields. Using (|6.9|) 
one finds rather complicated formula for the particle's velocity (see Appendix 0): 



p\pi - Vi)^'H^+p'Pk - v'Vk - nv\pi - Vi) ^ 

\pKpi-ViW + n'ip^-v^)iPi-n) ^^-^^ 



Now, inserting ( |8.1| ) into H = pkv'' + 7 we obtain 



H(q, p; fields) = V7^2 + p2 _ p2 ^ ^g 2) 



Let us observe that for a free particle, i.e. e = 0, H = m and Vk = 0, formula (|8.1|) gives 
relativistic relation between particle's velocity and momentum: 

Pk 

Vk 



\/ m'^ + p2 ' 

and formula ( |8.2| ) gives relativistic particle's energy: 



E{p) = Jm^ + p 



9 Poincare algebra 

Let us define the laboratory-frame Poincare generators. They are given by the Lorentz trans- 
formation of the rest-frame generators := {T-C,V^), TZk and S"^. Obviously, laboratory- 
frame four-momentum is given by ( |6.9| ) and ( |6.10| ): 

/ := l{n + v'Vk) , 
Pk ■= 'yvkH + 'y''^a\Vi . 

The static moment and the angular-momentum s"^ with respect to the particle's position 
q = (q^) are given as follows: 

Vk := ^^{a-'y,ni + ^eMmv'S"' + qkp\ 
:= 73(a-i)"^^5'= - ^e'^^^VkTZi + e'^^^quPi . 

With these definitions one can prove the following 
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THEOREM 7 Laboratory-frame generators p", pk, and form with respect to the 
reduced Poisson bracket { , }* the Poincare algebra. 

This already proves that our formulation is perfectly consistent with the Lorentz invariance 
of the theory. 

10 Particle in an external potential 

Suppose now that the particle moves in an external (generalized) potential U = [/(q, q, t). 
Then the Lagrango-Hamiltonian is given by: 



(n + a'^Uk - UmS"") - U . (10.1) 



To pass into Hamiltonian description in the particle's variables one has to apply the same 
procedure as in the Section ^. The primary constraints (f)^^^ are obviously given by ( |3.19|) 
but the secondary ones read 



da) 



dq 

Obviously, (j)""^ for a = 1,2] A; = l,2,3 are second-class constraints and the constraint algebra 
has the same form as ( |6.7| ) with [3k replaced by given by 

/3r' = /5fc + e-V(a^^)\Q. , (10.3) 

where 

^ dU ddU , 

is a vector of the generalized force in the laboratory frame. The complete Hamiltonian H 
on the constraint manifold reads: 

H = Pkv^ + r^n + u = ^{n + v^Vu) + u- v'^-^ . (10.5) 

Obviously, the constraint manifold V is endowed with a weak symplectic form Q and one 
can prove the following 

THEOREM 8 The Hamiltonian vector field Xu corresponding to l \10.d^ ) is well defined on 
the submanifold V* defined by the following non-homogeneous boundary condition 

DP(47rmvi>2 + e2vl/2,3)(0) = -e^^a^-^)''.Q,^ . 

r 

The proof is analogous to the proof of Theorem 
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Finally, the Hamiltonian structure for the particle in an external potential is a triple 
{V*, fl*, H), where Q* is a reduction of Q toV*. Due to the Theorem |]this structure is well 
defined, i.e. the initial data x; v) for the radiation field and for the particle uniquely 
determine the entire history of the system if the external potential is given. 

As an example consider the particle interacting with an external electromagnetic field 
/^^*. The generalized potential is given by: 

U{q, q, t) = eA-*(q, t) - eqA-*(q, t) , 

where Aq^^ and A*^^* stand for the four-potential of the external field in the laboratory frame. 
The generalized force ( |10.4| ) in terms of the laboratory-frame components Si and Bi of the 
external field now reads: 

Qi = e(Si{q,t)+e,jkv'l3\ci,t)) . 

In this case 

and the complete Hamiltonian reads 

H(q, v; fields) = 7 + v'^Vk) + eAg"* . 

Moreover, since U is linear in the velocity, we may easily express v in terms of q, p and 
fields (radiation and external). One obtains formula analogous to ( ^.1|) with pk replaced by 
Pk + eAl^^ which leads to the following expression for H: 



H(q, p; fields) = + (p + eA^-*)2 -V' + eA[;"* . 

Finally, let us observe that the "commutation relations" between particle's variables q, 
V and field variables (^I', x) have the same form as in Appendix |C| with f3 replaced by 
defined in ([Ta3|) . 



Appendixes 

A Canonical formalism for a 2-nd order Lagrangian 
theory 

Consider a theory described by the 2-nd order lagrangian L = L{q, q, q) (to simplify the 
notation we skip the index "z" corresponding to different degrees of freedom g*; extension 
of this approach to higher order Lagrangians is straightforward). Introducing auxiliary 
variables v = q we can treat our theory as a 1-st order one with lagrangian constraints 
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(p := q — V = on the space of lagrangian variables (g, q, v, v). Dynamics is generated by the 
following relation: 

dL{q,v,v) = —{pdq + n dv) = pdq + pdq + n dv + n di) . (A.l) 

where [p, tt) are momenta canonically conjugate to q and v respectively. Because L is defined 
only on the constraint submanifold, its derivative dL is not uniquely defined and has to be 
understood as a collection of all the covectors which are compatible with the derivative of 
the function along constraints. This means that the left hand side is defined up to fi{q — v), 
where /x are Lagrange multipliers corresponding to constraints = 0. We conclude that 
p = A is an arbitrary covector and ( |A.1| ) is equivalent to the system of dynamical equations: 






(-) 


d^ 


(-] 


dt 


[dv / 


dt^ 


\ dv 1 



dL dL dL 

dv dq dv 

The last equation implies the definition of the canonical momentum p: 

dL 

dv 

We conclude, that equation 

P 

is equivalent, indeed, to the Euler-Lagrange equation: 

5q dt^ [dv) dt[dv)^ dq~ ' ^ ■ ^ 

The Hamiltonian description is obtained from the Legendre transformation applied to ( [A.l|) : 

— dH = pdq ~ qdp + Tc dv — V dn , (A. 4) 

where H{q,p,v,n) = pv + tc v — L{q,v,v). In this formula we have to insert v = f,7r), 
calculated from equation tt = Let us observe that H is linear with respect to the 
momentum p. This is a characteristic feature of the 2-nd order theory. 

Euler-Lagrange equations ( [A.3|) are of 4-th order. The corresponding 4 Hamiltonian 



equations have, therefore, to describe the evolution of q and its derivatives up to third order. 
Due to Hamiltonian equations implied by relation ( [A.4|) , the information about succesive 
derivatives of q is carried by {v, 7r,p): 



V describes q 



dH 
dp 



hence, the constraint = is reproduced due to linearity of H with respect to p, 
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TT contains information about q: 



dH 



p contains information about Q 



dH dL 



TT 



dv dv 



the true dynamical equation equals 



P 



dH dL 



dq dq 



B Proof of Theorem 4 



Observe, that due to ( p.6| ) the field has the following expansion in the vicinity of the 
particle: 



An 



(B.i: 



where O(r^) denotes terms vanishing for r ^ like or faster. We shall prove that the 
Hamiltonian vector field Xh corresponding to the Hamiltonian ( |6.10| ) is well defined if and 
only if the fundamental equation ( |2.8| ) is satisfied. The field Xh is defined by 



XhJ n* 



6H 



iB.2) 



where f2* is given by ( |6.11D and Xh may be written in terms of coordinates as follows: 



^ dq^ ~'~ dv^ ~'~ ^ ~'~ ^ 6xa 



(B.3) 



Using (|B.2|) and (|B.3|) we get 

dH 



dq^ 



dH 

dv'- 



= Y' 

+ Za 

= -X 

+ Y'' 

- W 



d{lVk) da\ 



dv' 



dv' 



+ W 



6n 



SVn 



5n ^ „ 5Vn 
IVkT— + a fc- 



da\ 



d{-fVk) 
dv' " ' dv' 



(B.4) 



I n 



dv' 



dv^ 



Tin + I - 1 ^ 



dv' 



dv'' 



Za 



7a I- hU^ira 



^Xj 



(B.5) 
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^ n S'Pn 

6n 



+ a" 



7a k 



„ 6TZn , 5S 



6S' 



7a k-. + ^kra . 

dXA dXA 



+ Za 



Using the definition of H (see ( |6.10|) ) we liave 



^Xa 



0, 

l'^Vk{rL + v^Vi)+-iVk 
7 



7 



51-L , 5Vi 
h V — - 



Sxa Sx/ 



(B.6) 
(B.7) 

(B.8) 
(B.9) 
(B.IO) 

fB.in 



Now, calculating W"^ and Za from ( [B.7| ) and ( |B.6D respectively, inserting them to ( |B.5D and 
taking into account the Poincare algebra relations, we get 



-Ain + Bi^Vn-^^Si^iv'gnk[m 



Gvrrn 



X' 



CikH + Di^Vn - IsTia^Qnr m - 



67rro 



+ Y^[E,^^nn + FikmS"'] ,(B.12) 



where the following 3-dimensional objects depending upon the velocity v are introduced: 

Ai = -f\ - -f^aT^v'^gnk , 

^ 3„, „,ra ^ X ™ I ^ 2„n , , „,k^n 

k } 



Br 



Ik 



^ik 



Elk" 



Ikm 



-l{.9ki + I'^VkVi) + 7a" a\5(„^ , 
— dv^ + 7 ^'fca ; - cu^^a ^e^ 



'9(7a\) 9(7a" 



Irn \ 



Using the following properties of the function v'('r): 

2v9(r) - (1 -r)-i + V(r) 
2^'(r)-(l-r)-V(r)-(p^(r) 

and the identity 



v\eikiVyn + eumVk + eimkVi) = v'^ekir 
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one easily shows that 

Now, taking into account that 



7fi 



gik + 7 ^z^fc , 

we obtain finally the following equation for X'^: 

The matrix {g^i + •j'^VkVi) is nonsingular and its inverse equals to {g^'' — v^v^). Therefore 

Xk 2 / kl k l\ k 

= 7 ((? — V V )Vi = V . 

Now, let us consider equation ( [B.4|) and apply the same strategy as in the case of equation 
( |B.5|) . Calculate and Za from (|B.7] ) and ( |B.6|) respectively and insert them into 
Then, inserting X' = f ' and using once more the Poincare algebra structure one gets: 



\67rro 

l{9ki + 1^VkVi)H - a!'ia\gnr (u - m + 



(B.13) 



Now, observing that 
we finally obtain equation for F': 



m I = 7 '^a\ ( TT-— «n + ef3r 



V67rr, 







[B.14) 



The Hamiltonian vector field Xh is well defined if and only if the singular terms proportional 
to Tq ^ cansel out, i.e. 



eY\gki + I'^VkVi) = 7 ^a\a^ 



and therefore 



e e 



When the above equation holds then from ( p.l4j) 

- niY^gki + -f'^VkVi) = e-f-^a\P„ 

and finally 



which ends the proof of Theorem 0. 
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C Commutation relations 

The complete set of "commutation relations" read: 

{q\q'Y = 0, 

{v',v'r = 0, 

{g^xl(x)r = 0, 
{g^xI/2(x)r = 0, 

{^;^*^(x)}* = 0, 

Kxi(x)r = 0, 



and for field functional : 



{*^(x),*^(y)r = 0, 

{*'(x),xi(y)r = <^'(x-y), 
{xi(x),xi(y)}* = 0, 

{xi(x),x2(y)r = ^^Imm^^'^'^' 

{*^(x),*^(y)r = - f/rr n^(lyl) - /tt fT'^(W)^ 

e V<^X2(x) |y| 5x2(y) |x| / 

{X2(x),x2(y)r = f n^dyl) - 7^ n^(W)l 

2e V(^X2(x) |y| 5x2(y) |x| J 

{*'(x),x2(y)r = 
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Functionals fCi and A; are defined as follows: 



/C, := n^ + r^^v'eu^s'^ 



K := Vi + Vi[n + -K'l3i 

m 



e A , X, 



D Derivation of the formula v = v(p). 

From ( |6.9| ) we have 

Ck := Pk-Vk= {^n + ^ v'Vi) Vk . (D.l) 



Now, let X = v'^, A = Hy/ckC^ = 7i|c|, B = c{P^ and C = cipK From ( p. ID vectors c = (ck) 
and V are parallel and therefore 



Vk = -r-rCk ■ 
c 



Multiplying both sides of ( |6.9| ) by we obtain 

p,v'' = ^{n + Vkv'') . 

Moreover, 

Vkv" = Prc'Vk = . 
|c| |c| 

and 

k W\ k 
PkV = -Ac^Pk = ric . 

\c\ \c\ 

Therefore, inserting ( p.4|) and ( p.5|) into ( p.3| ) we obtain 

' 1 



C = A, 



X 



+ B- 



l — X a/1 — X 
The above equation is a square equation for ^/x■. 

{A^ + C^)x + 2AB^ + - C') = . 
The discriminant A for ( p.7| ) equals: 

A = 4^2(^2 + 52 - C^) . 
Calculating x and inserting into ( D.2 ) we finally obtain (|8.1|) . 



(D.2) 



(D.3) 



(D.4) 



(D.5) 



(D.6) 



(D.7) 
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